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It is a well-established fact that massive cosmological objects exhibit a "geometrical bias" that 
boosts their spatial correlations with respect to the underlying mass distribution. Although this 
geometrical bias is a simple function of mass, this is only half of the story. We show using numerical 
simulations that objects that are in the midst of accreting material also exhibit a "temporal bias," 
which further boosts their clustering far above geometrical bias levels. These results may help to 
resolve a discrepancy between spectroscopic and clustering mass estimates of Lyman Break Galaxies, 
a population of high-redshift galaxies that are caught in the act of forming large numbers of new 
stars. 



I. INTRODUCTION 

Large-scale structure in the Universe is believed to 
have originated from a primordial Gaussian random 
field of matter fluctuations, the product of quantum 
fluctuations that were shifted to larger spatial scales 
during cosmological inflation [e.g. 1]. Cosmic Mi- 
crowave Background observations show that when the 
Universe was 100,000 years old, the gaseous compo- 
nent was extremely smooth, with temperature vari- 
ations ^ 10~^. These tiny inhomogeneities, in con- 
cert with inhomogeneities in the unseen massive dark- 
matter, were amplified through gravitational insta- 
bility, eventually forming the galaxies, clusters, and 
other cosmological objects we see today. 

While the precise details of structure formation are 
highly complex, the simplicity of the initial random 
field allows us to easily compute the overall distribu- 
tions of cosmological objects. Galaxies and galaxy 
clusters represent the peaks in the initial density dis- 
tribution which accrete matter at the expense of the 
diffuse regions between them, and thus their number 
densities can be simply related to the number den- 
sities of the peaks of the initial random fleld. This 
technique has been applied most cleanly to galaxy 
clusters, whose densities and evolution provide strong 
constraints on the overall matter density 0|- 

Similarly, because the peaks in a random field are 
more clustered than the overall distribution, the spa- 
tial clustering of cosmological objects is stronger than 
the underlying mass distribution. Furthermore, this 
"geometrical bias" is a systematic function of the 
mass of these structures, an effect that has been well- 
studied analytically and numerically [e.g. 3, 4, 5]. 

Yet, this is only half of the story. Here we conduct 
a detailed numerical simulation that shows that the 
spatial correlation function of objects that are in the 
midst of accreting substantial amounts of material is 
significantly enhanced over that of the general popu- 
lation. This temporal bias causes them to mimic the 



properties of higher-mass structures, with important 
astrophysical implications as discussed below. 

The structure of this work is as follows: In §2 we 
describe our numerical simulation, discuss our group- 
finding algorithms, and develop a robust definition of 
accreting groups. In §3 we present our results for the 
spatial correlation functions of these samples, and in 
§4 we discuss the astrophysical implications of our re- 
sults. Further details of this study are given in 



II. SIMULATIONS AND GROUP FINDING 

Our numerical simulation traced the growth of pri- 
mordial density fluctuations by dynamically evolving 
a large number of point test particles. The distri- 
bution of these particles was then used to determine 
the nonlinear evolution of the spatial correlation func- 
tion of massive objects at late times. Driven by mea- 
surements of the Cosmic Microwave Background, the 
number abundance of galaxy clusters, and high red- 
shift supernova estimates [e.g. 7, 2, 8] we focused our 
attention on a Cold Dark Matter cosmological model 
with parameters H — 70 km/s/Mpc, = 0.3, f^A = 
0.65, Qb = 0.05, and erg = 0.87, where H is the Hubble 
constant today, Qq, fl^, and Qb are the total matter, 
vacuum, and baryonic densities in units of the critical 
density, and o"! is the present variance of linear fluc- 
tuations on the 8 x (100/70) Mpc scale (where 1 Mpc 
is 3.26 X 10^ fight years). 

Periodic boundary conditions, which approximate 
large-scale homogeneity and isotropy, were taken, and 
the mass within our simulation volume was held fixed. 
Thus the overall box expanded along with the cosmo- 
logical expansion, such that each side at any given red- 
shift z was 73/(1 -I- z) Mpc across. This box was pop- 
ulated with 350'^ dark matter particles that interact 
only gravitationally and represent the dominant mass 
component of the Universe. The mass of each particle, 
4.3 X 10^ solar masses (Mq), was chosen to match the 
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observed mass density of the Universe, and the simu- 
lation was started at an initial redshift of z = 49. The 
simulation used a parallel OpenMP-based version of 
the HYDRA code [9, 10] with 64-bit precision. 

To demonstrate the robustness of our results we 
have chosen two distinct group-finding approaches; 
the friends-of-friends approach (FOF) and the 
HOP algorithm FOF works by linking together 
all pairs of particles within a fixed "linking-length" 
of each other, and then taking each such group of 
"friends" to be an identified cosmological object. Al- 
though it remains popular, the FOF mass estimates 
are known to have significant scatter due to a problem 
that can occur as small strings of particles fall within 
the linking length. 

The HOP algorithm works by using the local den- 
sity for each particle to trace ('hop') along a path of 
increasing density to the nearest density maxima, at 
which point the particle is assigned to the group de- 
fined by that local density maximum. As this process 
assigns all particles to groups, a 'regrouping' stage is 
needed in which a merger criterion for groups above 
a threshold density Souter is applied. This criterion 
merges all groups for which the boundary density be- 
tween them exceeds S saddle, and all groups thus iden- 
tified must have one particle that exceeds Speak to be 
accepted as a group (see for explicit details). 

Beginning from z = 4.89, we saved particle posi- 
tions every 50 million years up to the final output at 
z = 3. For the final 5 outputs we found FOF groups 
using a linking parameter of 6 = 0.18, and HOP 
groups using the parameters: Ndens = 48, N^op — 20, 

Emerge = 5, Speak = 160, Swaddle = 140, and Souter = 

80. Visual inspection showed strong similarities be- 
tween the two populations, with a small amount of 
unavoidable noise coming from groups around the 80 
particle resolution limit (a group found by FOF at this 
limit may not be found by HOP and vice versa). We 
compare groups from one output to another by trac- 
ing back all particles with a given group index from 
the later output to the earlier output. Particles that 
show no membership to a group at the earlier time 
are regarded as 'smooth infall' while other non-null 
indices describe the merger history of the object. 

To give a rough estimate of the accuracy of the 
group finding methods in Fig. 1 we plot the mass of 
the most massive progenitor at ti{z = 3.059), versus 
the mass at t2{z — 3), such that At — h x 10^ years. 
As compared to the FOF groups, the smaller frac- 
tion of HOP groups lying above the equal mass line 
shows that the HOP algorithm identifies groups that 
are more likely to be massive at later outputs. The 
effect of this difference is significant. 

Our definition of accreting groups is similar to that 
of except that we select the subset that grew 
by 20% from t-^ to ^2, which implicitly includes mass 
accretion via smooth infall and results in 545(980) 
HOP (FOF) groups ff At = 5 x lO"^ years. Note that 
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FIG. 1; Comparison of group growth. The FOF algorithm 
exhibits a significant amount of scatter in mass estimates 
between outputs. Only 67% of groups grow from time t\ 
to t2, compared with 82% for HOP. 

the mass of each group is that at the end of each time 
interval, such that we tag all groups that experienced 
appreciable infall. The 20% value is arbitrary, but we 
selected it primarily because it appears to lie outside 
the central 'noise' band in the FOF data (see Fig 1). 
The groups corresponding to this cut appear in Fig 1, 
as points to the right of the dashed lines. 

III. TEMPORAL BIAS 

In Figure|21we show the spatial correlation functions 
of the groups selected by both the HOP and FOF al- 
gorithms and compare them with ^(r) of the accreting 
groups. This function measures the excess probability 
of finding a pair of groups at a given separation r rel- 
ative to a random distribution, and is calculated for a 
separation bin ri as 

i-KC(n) = iV(n)/iVra„dom(n), (i) 

where N{ri) is the number of pairs separated by 
distances between n and n+i, and A'random(?';) = 
^N'^^{rf^^-rf)/V, with N the total number of 
groups and V the volume of the simulation. In 
the accreting case we co-added the correlation func- 
tions calculated from the differences from the last 
four At = 5 X 10^ year intervals and the last two 
At = 10 X 10'' year intervals. Radial bins of 1/80 the 
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FIG. 2: Spatial Correlation Functions. In each panel the dashed line shows the correlation function for all the groups, 
while the points connected by the solid line show ^(r) for groups that have accreted appreciable msiss in the last At 
years. Panels are labeled by their mass range and At values, and in each panel, the dotted line shows the correlation 
function of all the groups in the next highest mass bin. The top two rows were generated from a set of groups selected by 
the HOP algorithm, while the groups in the lower two rows were selected using the FOF approach. The shaded region in 
the central panels represents the observed correlation function of TZab < 25.5 Lyman break galaxies as computed in fl^ 
by inversion of the angular correlation function. A 10% accretion threshold is applied in the IO^^'^A^q case to increase 
the number of measured groups. 



simulation size, corresponding to 0.92 comoving Mpc, 
were taken throughout. For comparison, in each panel 
of Fig. 12 we also show the correlation function of all 
the groups in the next largest mass bin. The ampli- 
tudes of the correlation functions obtained using the 
full set of HOP and FOF groups agree with each other 
to within statistical uncertainties, as well as with an- 
alytical estimates. 

The upper row demonstrates a clear enhancement 
of the clustering of accreting groups at both the 
10^°-^Mq and IO^^^Mq mass scales, with their cor- 
relation functions roughly matching those of objects 
three times greater in mass (no conclusion can be 
drawn from the high mass bin as the sample is too 
small). This "temporal biasing" arises from the fact 
that both objects accreting substructure as well as 
those experiencing considerable smooth infall tend to 
be found in the densest regions of space, which are 
themselves highly clustered. This conclusion is sup- 
ported by the fact that the average local overdensity 



of groups in the 10^°-^{10^^)Mq mass bin is 0.82(0.87) 
(measured in 4 Mpc (comoving) spheres, correspond- 
ing to a mass scale of 1.2 x 10^^ -M©), whereas the 
same mass bin for the entire population exhibits an 
overdensity of 0.60(0.73). 

In the second row of Fig. |21 we take a longer in- 
terval of At = 10 X 10^ yr. This has a only slight 
dampening effect on temporal bias, which can not be 
definitively distinguished from statistical noise in our 
measurements. In the At = 20 x 10^ yr case, however, 
only a very weak enhancement of ^(r) was measured. 

In the lower two rows of this figure, we repeat our 
analyses using the FOF group finder. Although this 
approach is more susceptible to statistical noise, the 
same trends are apparent as in the HOP case. If 
At = 5 X 10^ yr, this temporal bias is roughly equal 
to the geometrical bias of the groups three times more 
massive, while if At = 10 x 10^ yr, ^(r) is boosted to 
a slightly lesser degree. 

Finally, to quantify our results, we have computed 
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the effective temporal bias in each mass bin, At, 
and group finder. We define as the ratio of the 
correlation function of the accreting groups to the 
overall correlation function, weighted by the num- 
ber of points in each bin in the overall function; 

6j = X]?=o iT''"'' ' '^'^^'^^ carried out 

over all bins within r < 20 comoving Mpc. These val- 
ues are labeled in each panel, and in the Ai = 20 x 10^ 
yr case, 6^ ^ 1.1(1.0) in the 10^°-'%10^^-°)Mq HOP 
bins and 1.1(1.3) in the respective FOF bins. 

IV. ASTROPHYSICAL IMPLICATIONS 

While temporal bias is a general property of the 
peaks of a gravitationally amplified Gaussian random 
field, our results have specific implications for the 
large sample of z ~ 3 galaxies made available by the 
Lyman-break color-selection technique [T^ . Lyman 
break galaxies (LBGs) are observed to have enormous 
star- format ion rates on the order of ^ 50Mq per year 
[Tsl . implying that these objects are likely to be ac- 
creting large amounts of material. Furthermore, al- 
though the clustering of LBGs brighter than TZab < 
25.5 is roughly that expected from the geometrical 
bias of 10^^ Mq objects [e.g. 16, 17], the linewidths 
measured from a spectroscopic sample these galaxies 
correspond to total masses < IO^^Mq [l^ . 

To relate our result to LBGs we plot the spatial 
correlation function of TZab < 25.5 LBGs, as derived 
in 1^], in the center column of Fig. |21 Although 
there are significant uncertainties involved in comput- 
ing this quantity, since comparisons are more natu- 
rally conducted in angular coordinates, the shaded 
regions provide a guide to the range of ^(r) values 
consistent with observations. In these panels, we see 
that if At = 5 X 10^ yr is chosen, then temporal bias 
boosts the correlation function of IQ^^Mq groups into 



reasonable agreement with observations. 

This mass is marginally consistent with the upper 
mass bound inferred from the rotation curves of a 
somewhat bright {TZab 24) spectroscopic subset of 
LBGs Furthermore, only ~ 4% of all groups 

exhibit appreciable accretion in each At = 5 x 10^ 
year time interval and the density of IO^^Mq groups 
is ~ 2 X 10~^ Mpc'^, at 2; = 3 in our assumed cos- 
mology. Thus associating such objects with 5 x 10^ 
year starbursts results in a density 5 x 10"" Mpc^, 
comparable with that observed. 

While quite suggestive, these comparisons are not 
meant as a complete model, and may not prove to be 
the final explanation of the discrepant mass estimates 
of LBGs. Kinematic models have been explored, for 
example, in which the observed velocity dispersions of 
LBGs are much less than the circular velocities of the 
groups in which they are contained ■ What is clear 
however, is that this bias can not be ignored and must 
be carefully considered when interpreting the cluster- 
ing of these objects. While perhaps only part of the 
story, temporal biasing represents an important fac- 
tor that must be taken into account when studying 
the properties of Lyman break galaxies. 
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